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Abstract 

Predictions of the QLV (Quasi-Linear Viscoelastic) constitutive law are 
compared with those of the ABAQUS viscoelastic model for two simple 
motions in order to highlight, in particular, their very different dissipation 
rates and certain shortcomings of the ABAQUS model. 



1 Introduction 



As the demand for fuel efficient cars intensifies, tyre manufacturers try harder to 
estimate accurately the energy losses of their products. However the mechanics of 
tires are a most complex topic due to the variety of mechanical factors involved, 
such as nonlinear responses, structural inhomogeneities, complicated geometries 
and solicitations, hypo-, visco-, and hyperelasticity, etc. Consequently, advanced 
finite element codes are often called to the rescue in order to simulate real-world 
situations, and to evaluate the energy efficiency of a tyre. These commercial codes 
are often used as a "black box", and the validity of the answer is rarely put into 
question, even though it might provide a decisive argument in favor of, or against, 
the viability of a given tyre model. 

Here we examine the current implementation of nonlinear viscoelastic ef- 
fects in the Abaqus Finite Element Analysis (FEA) package. We find that it 
is not sound physically, and that it gives results that are not consistent with the 
thermomechanically-based Quasi Linear Viscoeleastic (QLV) model. We present 
both models in the next section, and highlight their main differences. Then we 
investigate two toy experiments, using the incompressible neo-Hookean material 
for the elastic response. We consider the equi-biaxial extension test in Section 
3, and the simple shear test in Section 4. In both cases, the Abaqus FEA model 
predicts higher energy losses than the QLV model. 



First we introduce some notations. We call F the deformation gradient; it is de- 
fined as F — dx/X, where x is the coordinate in the current configuration of a 
particle at X in the reference configuration. We call / its determinant: / = detF, 
and B = FF^ , C = F^F the associated Cauchy-Green strain tensors, where a 
superscript T denotes the transpose. We also introduce the velocity gradient D = 
[FF^^ + {FF^^Y]/2, where the superimposed dot represents the time derivative. 
When a solid is incompressible, every deformation is isochoric, so that 



We recall that the dissipated power per unit volume is G{t) •D{t) = tr {a{t)D{t)), 
so that E^, the dissipated energy per unit volume over a period T, is given by 



2 Abaqus FEA finite viscoelasticity model 



7=1, 



trD = 0. 



(1) 




(2) 
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Next we present the Abaqus FEA model. Section 4.8.2 of the Abaqus Theory 
Manual (Hibbit et al. 2007) gives the following constitutive relation to model 
nonlinear viscoelastic effects, 



a{t) = ae{t) + SYM{F{t 



F-\t] 



(3) 

where Cg is the instantaneous Cauchy stress response (elastic response at very 
short times) and G is the so-called memory kernel, characterizing the stress relax- 
ation (with G(0) = 1). Also, "SYM" denotes the symmetric part of the bracketed 
term; hence D = SYM{FF^^}. This constitutive relation is valid for compress- 
ible as well as incompressible solids, because in that latter case the hydrostatic 
term —pi in Ce (where p is a Lagrange multiplier) does not produce work, nei- 
ther in the instantaneous response, nor in the history term, as expected. 

For incompressible solids, 7 = 1 at all times and has the general form: 



(4) 



where i/Aj , 1//2 are scalar functions of t and of the first and second principal strain 
invariants. Then ([3]) reduces to 

a{t) = -p{t)I+\ifi{t)B{t) + \if2{t)B{tf 
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+ £ SYM |f(0 j'^ G{t - s) ^ri{s)C{syds 



F-\t)\, (5) 



where p{t) = p{t) + jQG{t — s)p{s)ds is arbitrary and remains to be determined 
from initial/boundary conditions. 

The Abaqus FEA model is reminiscent of, and similar to a well-established 
model of finite viscoelasticity, namely the Pipkin-Rogers model (Pipkin and Rogers 
1968), which reads in the incompressible case as (Wineman 1972, Johnson et al. 
1996) 

a{t) = -p{t)I + F{t)S^R[C{t),0]+ l^^^^{R[C{s),t-s])ds'jF{tf. (6) 

Here /> is a Lagrange multiplier resulting from the internal constraint of incom- 
pressibility and i? is a strain dependent tensorial relaxation function, with the gen- 
eral form 



R = (l)oI + (j)iC + (l)2C^, 



(7) 
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where ^2 are scalar functions of t and of the first and second principal 

strain invariants. 

By an appropriate choice of the the Pipkin-Rogers model reduces to the 
so-called Quasi-Linear Viscoelastic (QLV) model, which has proved to be a most 
successful phenomenological model for the behavior of non-linear viscoelastic 
solids, see references in (Johnson et al. 1996, Wineman and Rajagopal 2008). 
Equation ([6]) is derived rigorously by successive approximations from the basic 
physical requirements governing the behavior of solids with memory (such as the 
principle of determinism and local action, the principle of material objectivity, 
etc.), see the review (Drapaca et al. 2007). 

However we notice upon inspection of ([5]) and ([6]) that there are at least two 
differences between the models suggesting that the Abaqus FEA model does not 
rely on physical principles. First, the integral term in equation ([3]) is generally non- 
symmetric, in contrast with the integral term in equation ([6]). This is taken care of, 
in a somewhat arbitrary and ad hoc manner, by using the "SYM" operator. Also, 
the history (time integral) term in the Abaqus FEA model terminates with F{t)~^ 
in contrast with the history term in the the QLV model, which terminates with 

To emphasize the differences between each model, we henceforth focus on 
a viscoelastic incompressible solid which has an instantaneous response modeled 
by the neo-Hookean stress-strain relationship (already implemented into Abaqus), 

ae = -pI + ^oB, (8) 

where /io > is a constant, the initial shear modulus. Also, the time relaxation of 
the solid is assumed to be governed by a one-term Prony series expansion, so that 
after an infinite time the shear modulus settles at the value /ioo, say: 



— 4-1 1 _ ^ ) Q-'/'^ 



G{t) = ^+(l-^\e-'^\ (9) 
i"o V Mo y 

In that case, the Abaqus FEA model gives the identifications i/Zi = /io, V^2 = 0, 
so that 



a{t) = -p{t)I+^B{t) + {p^-^)SYM{F{t) 







C{s)ds 




f-Ut 



(10) 

In contrast, the QLV model gives the identifications ^0 = IM)G, ^1 = 0, ^2 = 0, 
and the history term turns out to be integrable, to give eventually 

a{t) = -p{t)I+[^^ + (Mo - ^o.)c-'/']B{t). (11) 
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Here the viscoelastic effect is at its most limpid, whereas it seems quite compli- 
cated in ([To]). We now look at two simple examples of motions. 



3 Equi-biaxial extension 

The equi-biaxial extension test of an incompressible solid is described by the fol- 
lowing motion, 

xi=X{t)Xu X2 = X{t)-^'^X2, X3 = A(0-i/2X3, (12) 

where X{t) is the stretch ratio in the direction of extension. Then clearly the 
deformation gradient is diagonal. 



F(0=Diag[A(0,A(0-^/^A(0-^/2 



(13) 



Assuming now that the solid is subject to uni-axial tension: aju] ^ 0, C7[22] = 
(7[33] = 0, allows us to compute the Lagrange multiplier p. Then we find the 
resulting Cauchy stress non-zero component is 

(14) 

for the Abaqus FEA model, and 

<y[n]{t) = [li^ + {lio-li^)e''/'] [X^it)-X{t)-'] , (15) 

for the QLV model. The difference between each model is now clearly apparent, 
and it reflects on the dissipated power, as we know confirm numerically. 

From an experimental point of view, it is common practice to have a dynamic 
displacement superimposed on a large static deformation. Here we consider that 
the neo-Hookean solid is strained in equi-biaxial tension by 30% from time t = 
to time t — I, and then made to oscillate with an amplitude of 20%: 

ri+0.3,, 0<,£1, 

[l.3 + 0.2sin£o(?-l), t>\. 

For the remaining parameters we pick 

lioo/lio = 0.5, T = 0.1s, (0=10.0s"\ (17) 
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Figure 1: Variations of the axial tension o^^^ with time, for the Abaqus FEA 
model (thick curve) and for the QLV model (thin curve) in the case of equi-biaxial 
tension. 

or we keep (O as a free parameter when we investigate the frequency dependence 
of the dissipated energy. 

Hence Figure [T] shows that at those values, the predictions from the Abaqus 
model differ widely from those of the QLV model. 

The numerical simulations are performed with Abaqus 6.7-1 on a bidimen- 
sional square element of 1 m side; we use a single CPS4 element (4 nodes, bilin- 
ear, plane stress) with an implicit solution scheme. Figure [2] shows clearly that the 
Abaqus model overestimates the energy dissipation compared to the QLV model. 



4 Simple Shear 

We now consider a simple shear of amount 7 say, 

Xi=Xi+ 7(0^2, X2 = X2, X3=X3, 

for which the deformation gradient is not symmetric. 



(18) 



F(t) 



I 7(0 
1 




(19) 



Simple calculations reveal that a sheared solid described by the Abaqus FEA 
model or by the QLV model is in a state of plane stress (aj,^] 7^ for /, j = 1,2; 
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Figure 2: Energy dissipation per unit volume over a single period for the Abaqus 
model (gray, solid curve) and for the QLV model (black, dashed curve) in the case 
of equi-biaxial tension. 

<7[3y] = for 7 = 1 , 2, 3) when the Lagrange multiplier p is taken as 

p{t) = jUoo - (jUoo - jUo)e"^/^. (20) 

Then we find for the Abaqus FEA model, 

g-(r-5)/T 



rt Q-(t-S)/T: 

= (Moo-mo) — - — r{s)r{t) ds, 



ty[i2](0 = Mor(0 + (Moo-juo)_^ ^ y{s) [2+r{s)r{t)-r{tf]ds, 



<y[22]{t) = (jUoo-jUo) 

and for the QLV model 



1 p-{t-s)/T: 





= [m» 


- - /io)e 


<7[12l(0 




- (jUoo - jUo)e 


<^[22](0 


= 0. 





-y{s) [y{s)-Y{t)]As, 



(21) 



(22) 

These latter expressions are intuitively expected for the neo-Hookean solid, be- 
cause for its instantaneous response at very short times, the Cauchy stress has 
components <7g[ii] = jUqT^, <7e[i2] = Mo7> <^e\22] = 0. In contrast, the Abaqus FEA 




Figure 3: Variations of C7[i2](0/Mo for the Abaqus analytical model (thick curve) 
and the QLV model (thin curve). 



model gives rise to a 0^22] component purely due to viscoelastic effects. This 
behavior is clearly non-physical, for in the case of a quasi-static loading, no trac- 
tion a[22] is required to shear a neo-Hookean solid, and yet the relaxation process 
would create such a component ex nihilo ! These discrepancies have grave conse- 
quences for the dissipated power because here 



D{t) 





7(0/2 




7(0/2 





so that (T-D = 7(?)a[i2](0, (23) 



which is obviously not the same depending on the model considered. 

We confirm these findings by testing the Abaqus FEA software against the 



current formulas (|2T])-(|23l). We take the amount of shear /(?) to vary as 



7(0 



0<t<l 
l + 0.2sin[(o{t-l] 



t > 1, 



(24) 



with the other parameters given by ( [17] ). Figures [3] and |4] display the variations of 



the 0^12] ^iid C7[j2] components computed from (|2Tj) and ( [22) ). Finally, we find that 
the Abaqus model overestimates the energy dissipation with respect to the QLV 
model, see Figure |5| Using the commercial software simulation, we recovered the 
thick curves, which confirms our opinion that equation ([3]) is actually implemented 
in the commercial code.. 
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Figure 4: Variations of C5'[22](0/Mo for the Abaqus FEA model (thick curve) and 
the QLV model (thin curve). The latter does not show clearly on the graph, be- 
cause it is equal to zero at all times. 
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Figure 5: Energy dissipation per unit volume over a single period in the case of 
the Abaqus model (gray, continuous curve) and of the QLVE model (black, dashed 
curve) for simple shear. 
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